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Wc present results for models of neutron stars and strange stars constructed using the 
Hartle-Thorne slow-rotation method with a wide range of equations of state, focusing 
on the values obtained for the angular momentum J and the quadrupole moment Q, 
when the gravitational mass M and the rotational frequency are specified. Building 
, on previous work, which showed surprising uniformity in the behaviour of the mo- 

ment of inertia for neutron-star models constructed with widely-different equations of 
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O . state, we find similar uniformity for the quadrupole moment. These two quantities, 

^ ' together with the mass, are fundamental for determining the vacuum space-time out- 

^ , side neutron stars. We study particularly the dimensionless combination of parameters 

QM/,P (using units for which c = G = 1). This quantity goes to 1 in the case of a 
Kerr-metric black hole and deviations away from 1 then characterize the difference 
between neutron-star and black-hole space-times. It is found that QMj for both 
, neutron stars and strange stars decreases with increasing mass, for a given equation of 

^SJ ' state, reaching a value of around 2 (or even less) for maximum-mass models, meaning 

0^ . that their external space-time is then rather well approximated by the Kerr metric. If 

' QM/J^ is plotted against compactness R/2M (where R is the radius), it is found that 

the relationship is nearly unique for neutron-star models, independent of the equation 
' of state, while it is significantly different for strange stars. This gives a new way of 

, possibly distinguishing between them. 
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1 INTRODUCTION 

Neutron stars are compact objects consisting of closely-packed neutrons together with protons and electrons (and at high 
densities also muons and possibly other particles such as hyperons, pions and kaons). Their central densities reach values where 
the microphysics is not well understood and so they could serve as laboratories for investigating the behaviour of nuclear 
matter under extreme conditions. A wide range of approaches has been used for modeling neutron stars and comparison 
between derived properties of the models and observations has enabled constraints to be placed on prescriptions for the 
nucleon-nucleon interactions (Haensel et al. 2007; Lattimer & Prakash 2007; Akmal et al. 1998; Rikovska Stone et al. 2003; 
Urbanec et al. 2010; Gandolfi et al. 2010; Newton & Stone 2009). An important, and widely discussed, alternative to the 
standard neutron-star picture is given by the possibility that some or all of the matter may consist of deconfined quarks. The 
most radical version of this, the strange star picture, envisages a star consisting entirely of deconfined quarks (apart, perhaps, 
from a thin crust of normal matter) and has its origin in the suggestion by Witten (1984) that such matter consisting of nearly 
equal numbers of up, down and strange quarks might represent the absolute ground state for strongly-interacting matter, 
even down to zero pressure. Strange-star models have been investigated and discussed by many authors (Haensel et al. (1986); 
Alcock et al. (1986); Alford et al. (1999); Farhi & Jaffe (1984); Colpi & Miller (1992) etc.). If these strange stars exist, they 
would have unusual properties because of being bound together by a combination of the strong force and gravity, rather than 
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just by gravity as is usually the case. It is then important to look for features of the internal and external space-time structure 
of compact stars, which might enable one to distinguish observationally between standard neutron stars and possible strange 
stars. 

Much of what we know about the rotation speeds of compact stars comes from observations of radio pulsars, thought 
to be rotating magnetized neutron stars emitting dipole radiation. Observed frequencies range from ~ O.lHz up to ~kHz 
(Manchester et al. 2005; Ghosh 2007). The fastest currently-known pulsar is PSR J1748-2446ad (Hessels et al. 2006) rotating 
with frequency of 716 Hz. Studying the rotational properties of compact stars is of great interest and several different 
approaches have been used for this within the context of general relativity. The first one involved the perturbative slow 
rotation approximation (taken to 2nd order in the neutron star angular velocity fi) as investigated by Hartle (1967); Hartle & 
Thorne (1968); Chandrasekhar & Miller (1974) and many others. The slow rotation approximation gives equations governing 
the structure of the rotating objects allowing one to calculate general properties such as the mass M, angular momentum 
J and quadrupole moment Q. These three parameters represent all that is necessary in order to characterize the external 
gravitational field within the slow-rotation approximation if one is retaining only terms in the perturbative expansion around 
the underlying non-rotating comparison object up to second order in the angular velocity Q.. Subsequently various finite- 
difference schemes have been developed for calculating models of rapidly rotating compact stars (Butterworth & Ipser (1975); 
Komatsu et al. (1989); Cook et al. (1994); Salgado et al. (1994); Stergioulas & Friedman (1995); Nozawa et al. (1998), etc.) 
and there is a freely downloadable code RNS (Stergioulas 1997) that can be used for calculating rapidly rotating neutron star 
models. Spectral methods have also been used and codes implementing these are available as a part of the LORENE package 
(Bonazzola et al. 1998). These have been used in various contexts (Lo & Lin (2011); Haensel et al. (2009); Bejger et al. (2010); 
Gondek-Rosihska et al. (2001); Amsterdamski et al. (2002), etc.) and have been compared with other numerical schemes by 
Stergioulas (2003); Berti et al. (2005); Nozawa et al. (1998). 

Properties of rotating compact stars were investigated in pioneering works by Hartle & Thorne (1968); Chandrasekhar & 
Miller (1974); Miller (1977) and others. Quadrupole moments were discussed by Miller (1977); Laarakkers & Poisson (1999) 
and others; recently Bradley & Fodor (2009) discussed the difference between the quadrupole moments of neutron stars and 
Kerr black holes. Some differences between rotating neutron stars and strange stars were pointed out by Bagchi (2010). The 
impact of the properties of neutron stars and strange stars on the behaviour of the external space-time and its astrophysically- 
important properties has been discussed by Torok et al. (2008); Haensel et al. (2009); Bejger et al. (2010); Gondek-Rosihska 
et al. (2001); Amsterdamski et al. (2002), etc. 

In previous works, several authors (Lattimer & Prakash 2001; Bejger & Haensel 2002) have shown that the factor I /MR^ , 
where / = J/Sl is the moment of inertia, follows an almost unique relation in terms of R/2M for most neutron-star equations 
of state whereas an entirely different behaviour was found for strange stars with an equation of state based on the MIT bag 
model. Motivated by these results, we present here a corresponding discussion for the quadrupole moments, finding again 
similar behaviour for the quantity QM/ . Throughout the paper, we use units for which c = G = 1. 



2 MODELS OF COMPACT STARS 
2.1 Equations of state 

We focus here on two families of compact objects. Standard neutron star models are taken to be composed of neutrons in 
/3-equilibrium with protons and electrons (and eventually muons and hyperons at higher densities) while strange star models, 
are taken to be composed of u,d and s quarks, (see, e.g. Haensel et al. (2007) for general overview of standard neutron-star 
models and Witten (1984); Haensel et al. (1986); Alcock et al. (1986) for details of the strange star hypothesis and models). 

We use here a wide range of equations of state for standard neutron-star matter, based on various different assumptions 
and methodologies. Variational theory is represented by the widely-used equation of state of Akmal et al. (1998): we use the 
model originally labelled as A18 -\- 5v -\- UIX* with the Argone 18 potential, including three-body forces and relativistic boost 
corrections. Relativistic mean field theory used to fit results obtained with direct Dirac-Brueckner-Hartree-Fock calculations 
is represented by the UBS equation of state (Urbanec et al. 2010). As representative of Brueckner-Hartree-Fock theory 
the equation of state of Baldo et al. (1997) labelled BBB2 is chosen. As a representative parametrization of the Skyrme 
potentials, we take the SLy4 equation of state - see e.g. Rikovska Stone et al. (2003) for detailed study of the Skyrme 
potential. As representative of the energy density functional developed by Bombaci (1995), model BPAL12 is used. We also 
include the relatively new equation of state of Gandolfi et al. (2010) based on the auxiliary field diffusion Monte Carlo 
(AFDMC) technique, which we label as Gandolfi. Equations of state including hyperons are represented by by those labelled 
as BALBNIHI (Balberg & Gal 1997) and GLENDH3 (Glendenning 1985)\ 

1 For models BALBNIHI, BBB2, BPAL12, BALBNIHI, GLENDNHl we used equation of state tables from the public domain code 
LORENE (Bonazzola et al. 1998). The APR equation of state was calculated using the effective Hamiltonian given in the original paper 



© 0000 RAS, MNRAS 000, 000-000 



Quadrupole moments of rotating neutron stars and strange stars 3 



As an equation of state for strange stars we are using the standard form of the basic MIT Bag model (Chodos et al. 1974) 
where the pressure P and baryon number density ub are related to the energy density £ via the equations 
1 



4(1 



4B), 
-2Qc/7r)^/" 



i£-B) 



3/4 



where B is the bag constant, whose value is related to the energy density of matter at zero pressure hy £ = 4B, and Qc is 
the strong interaction coupling constant. We here use the MIT Bag Model with standard parameters B = 10^*g.cm^^ and 
Qc = 0.15, and we also use twice this value for B in order to see its impact on the properties of the calculated models. The 
Bag Model gives a very simple representation of the matter in strange stars, but more sophisticated descriptions generally 
give very similar results for the main quantities of interest for this paper. 

We have chosen here equations of state representing different approaches to the microphysics of neutron star matter, but 
it is important to note that differences between different versions (or parametrizations) of the same approach can also be very 
significant, (see e.g. Rikovska Stone et al. (2003) for discussion of different parametrizations of Skyrme potential). We should 
stress here that our aim in selecting this set of equations of state was to choose a broad representative sample coming from a 
range of different microphysical approaches, without yet considering closer constraints coming from recent observations, such 
as those coming from the double pulsar (Podsiadlowski et al. 2005) , from the new highest-mass neutron stars (Demorest et al. 
2010; Antoniadis et al 2013) and from X-ray bursts (Ozel et al. 2010; Steiner et al. 2010). Our aim here is to show that even 
with taking this very broad set of equations of state, we see a surprising degree of convergence in some of the results obtained. 



2.2 The Hartle-Thorne method 

Hartle and Thorne (Hartle 1967; Hartle & Thorne 1968) developed a method for calculating models of rotating neutron 
stars using a slow rotation approximation that is valid for suitably small angular velocities with <^ GM/R^. It has been 
shown that this approximation can be used with good accuracy for almost all currently observed neutron stars, even for 
most millisecond pulsars (Stergioulas 2003; Berti et al. 2005). While it is, of course, of great interest to calculate very rapidly 
rotating models, the Hartle-Thorne method, when correctly used within its strict range of validity, is the most accurate 
method available, and it is appropriate for the vast majority of known neutron stars. 

The space-time metric around and inside rotating compact objects is given in the slow rotation approximation by the 
perturbation of a spherically symmetric metric retaining terms up to second order in the angular velocity Q (as measured 
from infinity); this can be expressed in the form 

f p2A 

ds" = -e^''[l + 2ho{r) + 2h2{r)P2]dt^ + e''Ul + —[2mo{r) + 2m2{r)P2] 

+ r^[l + 2k2ir)P2]{de^ + [d(P-u;{r)dtfsm'e}. (1) 

One can see that the perturbation away from the spherically symmetric non-rotating solution involves all of the metric functions 
fltti Srr, gee and g^^^ and that an additional term gt^ appears, describing the dragging of inertial frames. The perturbation 
functions ho{r), h2{r),mo{r),m2{r), k2{r) are quantities of order fl^ and are functions only of the radial coordinate r , while 
oj{r) is of order Q. The deviation away from spherical symmetry in the diagonal part of metric is given by the 2nd order 
Legendre polynomial 

^2 = P2(cos6i) = (3cos^6'- l)/2. (2) 

All of the perturbation functions must be calculated with appropriate boundary conditions at the centre and at the surface 
of the configuration. The second-order ones are labelled with a subscript indicating the multipole order of the perturbative 
quantities: I = for the monopole (spherical) deformations and / = 2 for quadrupole deformations representing the deviation 
away from spherical symmetry. Equations for the perturbation functions are derived from the Einstein field equations G^" = 
SttT^". If we put zero on the right hand side of these, we get the equations for the external vacuum space-time and then we 
can express the perturbation functions in terms of properties of the central object as measured by a distant observer: the 
mass M, angular momentum J and quadrupole moment Q (Hartle & Thorne 1968; Abramowicz et al. 2003; Chandrasekhar 
& Miller 1974). We take the energy-momentum tensor to be that of a perfect fluid 

T"" = (£• + P) [/" (7" -f Pg"" , (3) 



(Akmal ct al. 1998) and the UBS model represents the parametrization H in the original paper (Urbanec et al. 2010). The table for SLy4 
was kindly provided by Jifina Ri'kovska Stone. The Gandolfi equation of state was kindly provided by authors of original paper. 
^ The subscript in the metric functions u and A refers to the unperturbed Schwarzschild geometry. 
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where £ is the energy density, P is the pressure and [/'' is the fluid four-velocity. We deal with models rotating uniformly 
with angular velocity £7, and the non-zero components of the four-velocity are then 

V' = [-(ffii - 'i-'^gt.i. + '''' , u"* = w\ (4) 

The method used for calculating neutron-star models and obtaining the integral properties (M, J and Q), starting from 
specifying the central density and the rotation frequency, has been discussed by Hartle (1967); Hartle & Thorne (1968); 
Chandrasekhar & Miller (1974); Miller (1977). We use the same procedure here. 

We start by solving the equations for the unperturbed non-rotating objects. These are the equations of hydrostatic 
equilibrium, and the equations for the metric functions 

dP _ m(r) + 47rr^P 

dr " ' r(r - 2m(r)) ' 



dm 
dr 

du _ m{r) + 47rr^P 



47^£:r^ (6) 

(7) 



dr r(r — 2m(r)) ' 



e- = (^1-^j , (8) 

with appropriate boundary conditions at the centre of the star. For starting the integration, we need to have expressions for 
the values of the dependent variables a little away from the centre. These are obtained by making series expansions, leading 
to 

P Pc-(27r/3)(Pc+fc)(3Pc + fc)r', (9) 

m 4/3Tv£^r^, (10) 

u u, + {2n/3){Pc + £cy, (11) 

e^^ [1 ~ {8/3)TT£yy\ (12) 

The subscript c denotes the value of the quantity at the centre; Vc is initially unknown, but can be calculated by changing 
the variable to u — Vc and then matching the interior and exterior solutions at the surface of the star. 

For considering the rotational perturbations, we follow the earlier work in comparing rotating and non-rotating models 
having the same central density. Inserting the form of the Hartle-Thorne metric (1) and the energy momentum tensor for a 
perfect fluid (3) into the Einstein field equations then leads to differential equations for perturbation functions. We start with 
the one coming from the {t(f)) component 

(-«§)+ 4-^ 

where j = e^'^^"^' and uj — O, — lu. Near to r = 0, iD has following behaviour: 

Cj Cjc + ^ {£c + Pc) ^I'cr^ (14) 
5 

where (Zi^ is a constant whose value can again be found by matching with the exterior solution. Outside the star, one has 
2 / 

Cj{r) = n-^, (15) 

where the constant J is the angular momentum of the rotating neutron star (Hartle 1967). The angular momentum is then 
given by the formula 



2m(r) 



(16) 



T [d^J 

and both Cjc and J can then be calculated by matching Cj and dti/dr at the surface with their exterior solutions, if is 
specified. The moment of inertia / is given by the standard relation I = J/Q. 

The other field equations lead to the following differential equations for the / — and / — 2 perturbation functions. The 
/ = functions are given by 

— = {£ + P)—po + -r , (-J -3^..^, (17) 



dr dr 3 dr I r — 27Ti(r 



dpo dho _|_ 1 d /' r^pu^ 



mo(l + 87rr2p) ^■k{£ + Py^ 1 r^f ( dCj\'^ 1 d / r^fCj'^ , 



(r — 27Ti(r))2 r — 2m{r) 12 r — 2m(r) V dr / 3 dr Vr — 2m(r) 
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Near to r = 0, mo and po have the following behaviour: 



mo 



— (£ 

15 



• 2-2 5 



Po -^j'ci^cr^, 



while outside the star, where £ = P = 0, m{r) = m{R) = Mo and j = 1, they are given by 
mo 



SM 



(19) 
(20) 

(21) 
(22) 



r-2A/o r3(r - 2 Afo)' 

where 5Af is a constant giving the change in gravitational mass resulting from the rotation. For the I = 2 perturbation 
functions, where we use V2 — h2 + ^2 instead of ^2, we have 



dV2 



dr 

dh2 
dr 



dr 



dvo 



dr 



duj 



dr 



1 



3 - 2 d j 

" dT 



2U2 



r(r — 2m(r))dt'(j/dr 



-2 



di^o ^ 

dr 2(r — 2m(r))dfo/dr 



87r(£: + P) 



4m(r) 



dvo 



1 



dr 2(r — 2m(r))di^o/dr 
and m2(r) is given by 
m2 



Vdr7 3 



dfo , 1 

r 1 

dr 2(r — 2m(r))di'o/dr 



2 - 2 d j 
r w — — 
dr 



1 4 .2 (dCb\^ 1 3~2dj2 



3 ~2l 

r " dT 



r — 2m(r) 

The solutions of equations (23,24) are expressed as the sum of a particular integral and a complementary function 

h2 = 



V2 = 



where A is a constant to be determined and the homogeneous equations for complementary functions are 



dv. 



(C) 



dr 



dh. 



(C) 



„d//o, (C) 
^^dr''^ ' 

2v. 



(C) 



+ 



dr r(r — 2m(r))dz^o/dr 

Near to r = particular integrals have the behaviours 



2^ + 

dr 2(r — 2m(r))di'o/dr 



87r(£: + P) 



4m(r) 



,(C) 



'2 

, (P) 



br 



where a and b are constants which are related by 

b + Y(fc + 3P.)a = |^(£c + Pc)jl, 

while the complementary functions have the behaviours 



h. 



(C) 



Br' 



„(C) 



- — {£. + 2,P,)Br\ 



where B is another constant. The integrations are carried out with arbitrarily assigned values of a and B. 
Outside the star, h2 and V2 take the form 



h2 



V2 



J' 
J2 



\Mor 



+ K- 



1 

-* 

2Mo 



r 

Mo 



(23) 



(24) 



(25) 



(26) 
(27) 



(28) 
(29) 

(30) 
(31) 

(32) 

(33) 
(34) 

(35) 
(36) 



r- [r(r -2A/o)]i/2 

where K is & constant and the Qj^ are associated Legendre functions of the second kind (see equations (137) and (141) of the 
original paper Hartle (1967) for explicit formulae). The constant K can then be calculated by matching the internal values 
for 112 and /i2 to the external ones at the surface. 

Once these integrations have been carried out, the total mass of the rotating star is given as 
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Figure 1. Mass versus radius (Left) and compactness x = i?o/2Mo (Right) for non-rotating models with the selected equations of state. 



M = Mo + SM = Mo + mo(i?) + J7i?^ 
and its quadrupole moment is 



(37) 



(38) 



3 RESULTS 

3.1 Static non-rotating configurations 

We here denote integral quantities for the non-rotating object using the subscript (i.e. Mo is the mass of the non-rotating 
star and Ro is its radius) and we define tlie compactness parameter x of the star as its actual (circumferential) radius divided 
by its Schwarzschild radius 

The calculated mass-radius relations for the selected equations of state, for both neutron stars and strange stars, are 
plotted in Fig. 1 (Left) while Fig. l(Right) shows the resulting dependencies of mass on the compactness x. The vertical line in 
Fig. 1 (Right) corresponds to the condition Ro = 3Afo which is boundary of the situations where the neutron star has or does 
not have unstable circular null geodesies external to its surface. We can see that none of the selected equations of state allows 
the existence of "extremely compact objects" with Ro < 3Mo which would contain trapped null geodesies in their interior 
and thus could cool via the non-standard scenario suggested by Stuchh'k et al. (2009); however, the Gandolfi et al. (2010) 
equation of state allows configurations very close to that condition. The horizontal line represents the canonical neutron star 
mass M — 1.4 A/q. We can see that, excluding MIT2, this mass corresponds here to the interval « 2.4 to « 3.5 in compactness 
and so to a range of radii Ro of (4.8-7.0) Ado- Note also that according to Fig. l(Right), assuming the most astrophysically 
interesting masses of neutron stars, (M greater than around 1.25Mo), the corresponding compactness is x < 4 i.e. Ro < 8Mo. 

Concerning Fig. l(Left), we stress again that the range of equations of state taken here is extremely wide, and that some 
of these mass/radius curves are not consistent with the observational constraints mentioned earlier, although it should be 
pointed out that those constraints apply to certain particular objects and the equations of state which they seem to exclude 
could still be valid for other ones. It is by no means certain that one single equation of state (or microphysical model) would 
apply for all of these objects. 



3.2 Moments of inertia 

The moment of inertia of a neutron star model is calculated from the angular momentum using the standard formula I — 
Since the angular momentum J, as calculated here, is of first order in the angular velocity SI (with the rotational shape 
correction at order Q,'^ being neglected), dividing it by Q, leads to a quantity which does not vary with Q and depends only 
on the structure of the spherical non-rotating comparison object on which the perturbative expansion is based. It was shown 
by Lattimer & Prakash (2001); Bejger & Haensel (2002) that, for standard neutron stars, the dimensionless factor I/MoRo 



© 0000 RAS, MNRAS 000, 000-000 



Quadrupole moments of rotating neutron stars and strange stars 7 



0.6 



0.5 



0.4 



M o 



0.3 



0.2 



0.1 



Skyrme — 
UBS 
APR 
BBB2 - - 
BPAL12 - - 
BALBNIHI ■■ ■ 
GLENDNH3 - 
Gandolfi — 
MIT 
MIT2 — 
■ ■ ■ ■ I ■ ■ ■ ■ I 



1 



1.5 



2.5 3 3.5 

i?o/2Mo 



4.5 



0.25 



0.2 



0.15 



0.1 



0.05 



i i 1 i 1 1 1 1 i 


1 i i 1 1 i 1 1 i 1 i 1 i i 1 1 i 


i 1 1 i 1 1 1 1 i 1 i i 


APR 






UBS 






BBB2 - 







BPAL12 - 






- BALBNIHI 






iVlii 






MIT2 


















. . , . 1 , , , . 


. ... \ .... \ .... \ . . 


.. \ .... \ ... . 



1.5 



2.5 3 3.5 
i?o/2Mo 



4.5 



Figure 2. Left: The moment of inertia factor I /MqRq plotted versus compactness for tlie selected equations of state. The curve given 
by the "universal" formula is shown with the heavy solid line. Right: The angular momentum parameter j = J/M^ for objects rotating 
at 300 Hz plotted versus compactness x = Ro/2Mq. 



could be related to the compactness by a "universal" formula which is almost independent of the equation of state. In terms 
of our compactness parameter x — Ro/2Mo, 

= 029^' ^^3.33, 

I/MoRl = 2x + 3^^ ^< 3.33. (40) 
9a; 

This relation is plotted with the heavy solid line in the left frame of Fig. 2 and it can be seen that it fits rather well with the 
curves for the standard neutron star equations of state, despite the diversity of their formulations, whereas the curves for the 
strange stars are quite different (the two strange-star curves lie almost exactly on top of each other, despite having different 
values for the bag constant). The strange star moment of inertia factor I/MqRq for large values of x (i.e. for low mass strange 
stars), tends towards the value 2/5 which is the well-known result in classical physics for uniform density spheres. This is not 
surprising since these models do indeed have almost constant density profiles with the central energy density being close to 
the surface value, which also results in them having masses roughly following AIq oc Rq. 

The right-hand frame of Fig. 2 shows the angular momentum parameter j = J/Mq ( J/M^ taken to lowest order) plotted 
as a function of compactness for objects rotating at 300 Hz, which is around the maximum frequency for which the slow- 
rotation approximation can be considered as reliable for a neutron star with canonical mass and radius. (Note that this usage 
of j is not to be confused with that in section 2.2; we are caught here between two conflicting notations.) Since this j depends 
linearly on the rotational frequency, it is simple to calculate its value for any other rotational frequency, e.g. one tenth of the 
frequency would lead to one tenth of the value of j for a given equation of state and compactness^. 



3.3 Quadrupole moments 

The leading order contribution to the quadrupole moment (and the only one retained within our slow-rotation approximation) 
is of 2nd order in the angular frequency SI. A suitable dimensionless and frequency-independent quantity characterizing this 
is provided by the Kerr factor q — QMq/J^ {QM/ taken to lowest order) (Miller 1977) that can be used to represent the 
deviations of the external Hartle-Thorne metric away from the Kerr black hole metric (for which q = 1). For a given equation 
of state, the value of q is fully determined by the central density or pressure of the star or by its corresponding compactness 
X. Computed values of q for different degrees of compactness are shown in Fig. 3 for both neutron stars and strange stars 
and it can be seen that, once again, there is an almost universal relationship between q and compactness for the neutron star 
models, while the relation for the strange stars is quite different. 

We have again looked for a suitable analytic fitting formula to represent this almost universal behaviour for the neutron- 
star models. For doing this, we selected a linear dependency for the higher values of x and a quadratic one for the lower 
values, with the minimum corresponding to a- = 1, q = 1, i.e. we use the relations 

q — aix + ao, x > Xo, 



This statement is of course true if the compactness is replaced with any parameter of the non-rotating neutron star, e.g. with its mass 
Mo. 
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Figure 3. The Kerr factor q = QMq/J^ plotted against compactness for tlie selected equations of state. The approximate analytic 
relation is labelled as "fit" and is shown using the bold solid lino. 



q = b{x - 1)" + 1, X 4: xo, (41) 

where ai and ao are fitted parameters while b and xq are calculated assuming that the function is continuous and smooth at 
the point xo where the functions are matched. Under these circumstances, the relation for b is 

6 = aiV[4.(~ai-a() + l)l, (42) 
and the matching point xq is given by 

.0 = '-^^ - 1. (43) 
ai 

We have found that for ai = 3.64 and ao = —5.3, the analytic relation fits the calculated values very well, as shown in Fig. 3. 

It can be seen from Fig. 3 that q is systematically decreasing for increasingly compact models and seems to be tending 
towards the Kerr value g = 1 as the non-rotating comparison star gets closer to becoming a Schwarzschild black hole, i.e. as 
a:: — >■ 1 or i?o — !> 2Mo. This also corresponds to the stellar models getting closer to the maximum mass limit. The upper limit 
for q, for the most astrophysically interesting neutron-star models, is about 9. For the strange stars of given compactness, q 
is always larger than it is for neutron stars of the same compactness, but the tendency towards q = 1 as _Ro ^ 2Mo is seen 
for both families of objects. We find that q is almost identical for both values of the bag constant, similar to the situation for 
the moment of inertia factor. 

A key feature of the results is the systematic decrease of q as the mass increases and approaches its maximal value for 
any given equation of state, as shown in Fig. 4. We also note that for compact stars with the canonical mass Af = 1.4Mq, the 
value of q depends strongly on the assumed equation of state (as a result of the varying compactness of the models) and its 
value can range between q ~ 2 and g ~ 7 — 8. 
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Figure 4. The Kerr factor q = QMq j is plotted against the mass of neutron stars (and strange stars) for the selected equations of 
state. 



4 CONCLUSIONS 

In this paper, we have calculated models of rotating neutron stars and strange stars within the framework of the Hartle- 
Thorne slow-rotation approximation for a variety of equations of state, in order to investigate the behaviour of the parameters 
determining the external space-time. We have shown that for the most astrophysically interesting cases of objects with masses 
greater than ~ 1.25M0 , the compactness parameter x = Ro/2Mo is less than ~ 4 in the case of neutron stars and less than ~ 3 
for strange stars. For these models the angular momentum parameter j = J /Mq corresponding to the maximum currently- 
observed rotation speed of 716 Hz would come in the interval 0.2-0.5, if calculated with the slow-rotation approximation, 
depending on the mass and equation of state. For the lower part of this range, the slow-rotation approximation could be used 
consistently even for this rather rapidly rotating object. Most observed neutron stars can be accurately treated within the 
more-or-less analytic slow-rotation approximation without needing a more elaborate numerical treatment. 

Even if different equations of state for standard neutron-star matter give rather different neutron-star properties, some 
combinations of neutron-star parameters can be very accurately approximated using just the compactness of the neutron star. 
That this is so for the moment of inertia factor I /ME? was previously known; we have here demonstrated that the same is 
true for the Kerr factor q = QM / for which we have presented a new analytic fit. The interval for possible values of q for the 
most astrophysically interesting neutron-star models ranges from q ~ 1.5 for the most extreme objects close to maximal mass 
up to g ~ 9 for low mass objects. The analytical representations of the key parameters enable us to express the space-time 
metric around a slowly-rotating neutron star in terms of its mass M, radius R and rotational frequency /. 

The "universal" relation between the Kerr factor q and the compactness x of the neutron stars is almost independent of 
the equation of state. However, the equivalent relation between q and x for the strange stars is significantly different from this. 
If it becomes possible to constrain both q and x independently from observations, this could be used as a way of indicating 
whether strange stars may actually exist. We believe that this is another significant new result, to set alongside realising that 
for compact stars (both neutron stars and strange stars) with masses close to the maximum allowed for the given equation 
of state, the external space-time is very close to the standard Kerr space-time, a fact that can considerably simplify the 
modelling of accretion and optical phenomena in the vicinity of these compact stars. 
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The quadrupole moment can also play a key role in determining the position of the ISCO, especially for low-mass compact 
stars. This can have an impact on the validity of QPO models as discussed by Urbanec et al. (2010) and so it is necessary 
to include the quadrupole moment in the analysis of observational data, especially for masses that are close to the observed 
values (Torok et al. 2010, 2012). In a forthcoming paper, we plan to study the astrophysical consequences of our results and 
to compare them with those obtained by other methods. 

ACKNOWLEDGMENTS 

This work has been supported by the Czech grants MSM 4781305903, LC 06014 and GACR P209/12/P740. The au- 
thors further acknowledge an internal student grant from the Silesian University in Opava (SGS/1/2012) and the project 
CZ. 1.07/2.3.00/20.0071 "Synergy" in the frame of Education for Competitiveness Czech Operational Programme supporting 
the international collaborations of the Opava Institute of Physics. The work has also been supported by the CompStar research 
networking programme of the European Science Foundation. 

REFERENCES 

Abramowicz M. A., Almergren G. J. E., Kluzniak W., Thampan A. V., 2003, arXiv:gr-qc/0312070 
Akmal A., Pandharipande V. R., RavenhaU D. G., 1998, Phys. Rev. C, 58, 1804 
Alcock C, Farhi E., Olinto A., 1986, ApJ, 310, 261 

Alford M., Rajagopal K., Wilczek P., 1999, Nuclear Physics B, 537, 443 

Amsterdamski P., Buhk T., Gondek-Rosihska D., Kluzniak W., 2002, A & A, 381, L21 

Antoniadis et al 2013, to be published 

Bagchi M., 2010, New Astronomy, 15, 126 

Balberg S., Gal A., 1997, Nuclear Physics A, 625, 435 

Baldo M., Bombaci I., Burgio G. P., 1997, A & A, 328, 274 

Bejger M., Haensel P., 2002, A & A, 396, 917 

Bejger M., Zdunik J. L., Haensel P., 2010, A & A, 520, A16+ 

Berti E., White F., Maniopoulou A., Bruni M., 2005, MNRAS, 358, 923 

Bombaci I., 1995, in I. Bombaci, A. Bonaccorso, A. Fabrocini, et al. ed.. Perspectives on Theoretical Nuclear Physics An 

equation of state for asymmetric nuclear matter and the structure of neutron stars, pp 223-237 
Bonazzola S., Gourgoulhon E., Marck J., 1998, Phys. Rev. D, 58, 104020 
Bradley M., Fodor G., 2009, Phys. Rev. D, 79, 044018 
Butterworth E. M., Ipser J. R., 1975, ApJ Letters, 200, L103 
Chandrasekhar S., Miller J., 1974, MNRAS, 167, 63 

Chodos A., Jaffe R. L., Johnson K., Thorn C. B., 1974, Phys. Rev. D, 10, 2599 
Colpi M., Miller J. C, 1992, ApJ, 388, 513 

Cook G. B., Shapiro S. L., Teukolsky S. A., 1994, ApJ, 424, 823 

Demorest P. B., Pennucci T., Ransom S. M., Roberts M. S. E., Hessels J. W. T., 2010, Nature, 467, 1081 
Farhi E., Jaffe R. L., 1984, Phys. Rev. D, 30, 2379 

Gandolfi S., lUarionov A. Y., Fantoni S., Miller J. C, Pederiva F., Schmidt K. E., 2010, MNRAS, 404, L35 
Ghosh P., 2007, Rotation and Accretion Powered Pulsars. World Scientific, Singapore 
Glendenning N. K., 1985, ApJ, 293, 470 

Gondek-Rosihska D., Stergioulas N., Bulik T., Kluzniak W., Gourgoulhon E., 2001, A & A, 380, 190 

Haensel P., Potekhin A. Y., Yakovlev D. G., 2007, Neutron Stars 1: Equation of State and Structure. Vol. 326 of Astrophysics 

and Space Science Library, SPRINGER 
Haensel P., Zdunik J. L., Bejger M., Lattimer J. M., 2009, A & A, 502, 605 
Haensel P., Zdunik J. L., Schaefer R., 1986, A & A, 160, 121 
Hartle J. B., 1967, ApJ, 150, 1005 
Hartle J. B., Thorne K., 1968, ApJ, 153, 807 

Hessels J. W. T., Ransom S. M., Stairs I. H., Freire P. C. C, Kaspi V. M., Camilo P., 2006, Science, 311, 1901 

Komatsu H., Eriguchi Y., Hachisu I., 1989, MNRAS, 237, 355 

Laarakkers W. G., Poisson E., 1999, ApJ, 512, 282 

Lattimer J., Prakash M., 2007, Phys. Rep., 442, 109 

Lattimer J. M., Prakash M., 2001, ApJ, 550, 426 

Lo K., Lin L., 2011, ApJ, 728, 12 

Manchester R. N., Hobbs G. B., Teoh A., Hobbs M., 2005, Astron. J., 129, 1993 



© 0000 RAS, MNRAS 000, 000-000 



Quadrupole moments of rotating neutron stars and strange stars 



Miller J. C, 1977, MNRAS, 179, 483 

Newton W. G., Stone J. R., 2009, Phys. Rev. C, 79, 055801 

Nozawa T., Stergioulas N., Gourgoulhon E., Eriguchi Y., 1998, Astron. Astrophys. SuppL, 132, 431 
Ozel F., Baym G., Giiver T., 2010, Phys. Rev. D, 82, 101301 

Podsiadlowski P., Dewi J. D. M., Lesaffre P., Miller J. C, Newton W. G., Stone J. R., 2005, MNRAS, 361, 1243 
Rikovska Stone J., Miller J. C, Koncewicz R., Stevenson P. D., Strayer M. R., 2003, Phys. Rev. C, 68, 034324 
Salgado M., Bonazzola S., Gourgoulhon E., Haensel P., 1994, Astron. Astrophys. SuppL, 108, 455 
Steiner A. W., Lattimer J. M., Brown E. P., 2010, ApJ, 722, 33 

Stergioulas N., 1997, http://www.gravity.phys.uwm.edu/rns. UWM Centre for Gravitation and Cosmology 
Stergioulas N., 2003, Living Reviews in Relativity, 6 
Stergioulas N., Friedman J., 1995, ApJ., 444, 306 

Stuchh'k Z., Torok G., Hledi'k S., Urbanec M., 2009, Classical and Quantum Gravity, 26, 035003 

Torok G., Bakala P., Stuchlik Z., Cech P., 2008, Acta Astronomica, 58, 1 

Torok G., Bakala P., Sramkova E., Stuchh'k Z., Urbanec M., 2010, ApJ, 714, 748 

Torok G., Bakala P., Sramkova E., Stuchh'k Z., Urbanec M., Goluchova K., 2012, ApJ, 760, 138 

Urbanec M., Betak E., Stuchh'k Z., 2010, Acta Astronomica, 60, 149 

Urbanec M., Torok G., Sramkova E., Cech P., Stuchh'k Z., Bakala P., 2010, A & A, 522, A72+ 
Witten E., 1984, Phys. Rev. D, 30, 272 



© 0000 RAS, MNRAS 000, 000-000 



